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Abstract 

The derivation procedure of exact ground-states for the periodic Anderson 
model (PAM) in restricted regions of the parameter space and D = 2 dimen- 
sions using plaquette operators is presented in detail. Using this procedure, 
we are reporting for the first time exact ground-states for PAM in 2D and fi- 
nite value of the interaction, whose presence do not require the next to nearest 
neighbor extension terms in the Hamiltonian. In order to do this, a completely 
new type of plaquette operator is introduced for PAM, based on which a new 
localized phase is deduced whose physical properties are analyzed in detail. 
The obtained results provide exact theoretical data which can be used for 
the understanding of system properties leading to metal-insulator transitions, 
strongly debated in recent publications in the frame of PAM. In the described 
case, the lost of the localization character is connected to the break-down of 
the long-range density-density correlations rather than Kondo physics. 
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Typeset using REVT^ 



I. INTRODUCTION 



The periodic Anderson model (PAM) is one of the basic models describing strongly cor- 
related systems whose characteristics fit to be interpreted in a two-band picture^. The model 
is largely applied in the study of heavy-fermion systems^, intermediate- valence compounds^, 
or even high critical temperature superconductors^. In contrast however with other mod- 
els used in the understanding of phenomena created by strong correlation effects, where at 
least in one dimension exact solutions exist (for example the Hubbard modeP), the physics 
described by PAM is almost exclusively interpreted based on approximations. This is due 
to the fact that only few results are exactly known about the behavior of PAM. 

Indeed, what we know about the physical behavior of PAM in rigorous terms can be 
summarized as follows: The first results, related to the ground-state of decorated hyper- 
cubic lattices in the limit of infinite interaction strength has been provided by Brandt and 
Giesekus^ followed by a non-magnetic ground-state restricted to ID and on-site repulsion 
U — oo case obtained by Strack^. This solution has been extended to D — 2,3 as well, atU — 
qq8,9 Again for infinite on-site Hubbard repulsion has been demonstrated that at quarter 
filling the ground-state is unique with a defined total spin^°. It has been also underlined 
that the model becomes solvable in the case of constant and infinite range hopping^^. We 
further know, that for only on-site hybridization and without direct hopping in the correlated 
band: the symmetric half-filled case is spin-singlet^^ and in ID also pseudo-spin singlet at 
half filling anti- ferromagnetic correlations are present^^, and in ID, 2D long-range order of 
ferro, anti-ferro and pairing type is absent^^. Recently have been published the first exact 
ground-states at finite U in ID^^'^'' and 2D^^, respectively. Concerning 2D, the reported 
ground-states-*^^ require next to nearest-neighbor (NNN) one-particle terms as well in the 
Hamiltonian (H), and from the obtained solutions, especially the physical properties of the 
itinerant one has been described in detail. 

The deduction of exact ground-states in D = 2 dimensions is of large interest in general, 
and for strongly correlated systems in special. In the case of PAM, at finite and nonzero 
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value of the interaction, the only one procedure doable working at the moment in this re- 
spect, is based on plaquette operators introduced in Ref.^^, and a decomposition of the 
on-site Hubbard interaction as described in Ref.^^. During this procedure, H is such trans- 
formed to contain in a positive semidefinite form products of plaquette operators. Since 
in 2D the product of plaquette operators generates NNN one-particle terms as well, the 
general impression created by the method suggests that the applicability of the procedure 
is intimately connected to NNN contributions in 2D, and the deduced exact ground-states 
are fingerprints of this fact. 

In this paper, presenting for the first time exact ground-states for PAM in 2D at finite 
value of the interaction, in restricted regions of the parameter space and without NNN 
type of extension terms in if, we demostrate that the plaquette operator procedure and 
the ground-states deduced with it are not necessarily connected to the presence of NNN 
extensions in H. In order to clarify these aspects (i) the plaquette operator technique is 
analysed in detail in general terms and 2D, (ii) a completely new type of plaquette operator 
is introduced which allows the deduction of the presented results, (iii) the obtained new 
localized exact ground-state is analyzed and described in detail, and (iv) implications of the 
results relating the metal-insulator transition in frame of PAM are presented. 

The deduced new ground-state is a completely localized state. In order to characterize 
this phase, after obtaining the exact ground-state, all relevant ground-state expectation 
values and correlation functions have been exactly calculated and analyzed. The obtained 
state is paramagnetic, and based on a coherent control which it has on the occupation 
number of all lattice sites, it introduces long-range density-density correlations into the 
system, producing a localized state. 

Concerning implications to physical systems, we mention the intense activity in the field 
related to the understanding of the metal-insulator transition (MIT) in frame of PAM. The 
subject has an almost 30 years of history^^, and gained renewed interest in the last pe- 
riod based on the observed MIT similarities between the Hubbard model and PAM, used 
for example in the understanding of the iso-structural electronically driven MIT transitions 
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(like the 7 — > o; transition in Ce compounds) Since the exactly deduced ground-state 
energy values presented in this paper are not containing exponential factors characteristique 
to Kondo type of behavior, the results reported here underline that at least in some re- 
gions of the parameter space, a localization-delocalization transition in frame of PAM is not 
necessarily connected to Kondo physics. 

The remaining part of the paper is structurated as follows: Section II. presents the 
Hamiltonian we use. Section III. describes the plaquette operator technique. Section IV. 
presents the transformation of the starting Hamiltonian into a new expression containing 
plaquette operators. Section V. describes the detected new exact ground-states. Section VI. 
concludes the paper, and Appendices A.- D. containing the mathematical details of the 
starting points of the paper close the presentation. 

II. THE EXPRESSION OF THE HAMILTONIAN. 

We start with a generic PAM Hamiltonian written for 2D square lattice as 

H^fd + ff + Ef + V + U , (1) 

where, the contributing terms in order, are representing the kinetic energy of d electrons {T(i) , 
the kinetic energy of / electrons (T)), the on-site / electron energy (-B/), the hybridization 
{V), and the on-site Hubbard interaction written for / electrons (j — UUf, the last con- 
tribution representing the interaction term, and t/ > being considered during this paper. 
The presence of Tf in H is motivated by the overwhelming evidence, that the heavy-fermion 
materials contain in their real band structure around the Fermi level {Ep) very narrow, 
hybridized bands, which exist already at temperatures far above the thermodynamically 
determined Kondo temperature, being relatively T independent and holding an accentuate 
/ character^^. 

The interaction term during this paper is exactly transformed in the form 

Uf = T. MAi + E(E - 1) , (2) 

i i o- 
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where, the positive semidefinite operator P' — — — fi^| + n^|r7,^l^) defined by Eq.(2) 
requires for its lowest zero eigenvalue at least one / electron on every lattice site-*^^. As will 
be clarified in Section V., the representation presented in Eq.(2) is a key feature from the 
point of view of the interaction term in the process of the deduction of exact ground-states 
in the frame presented here. 

The hybridization V is considered to be build up from a local Vq, and a nonlocal Vni 
contribution, V — Vq -\-Vni. Thus, the local one-particle terms of the Hamiltonian are Ef 
and Vq whose expressions become 

Ef^E^Y. <a , Vo = Y.WlAa + H.c.) . (3) 

The non-local one-particle contributions remain to be presented. In order to make the 
notations clear, instead of the site-numbering notation (i) we use here the i vectorial notation 
for the lattice sites. The kinetic energy contribution f — fd + ff thus is given by 

= E ( dlA-,r,a + 4 fU+r,a + H.C. ) , (4) 
i,r,a 

and the non-local hybridization becomes 

Vnl^Ei Vf dljl+r,. + Vj" flJ,+r,. + H.C. ) . (5) 

The notation of the non-local hybridization matrix elements by the superscripts {df) and 
(fd) is given by mathematical convenience, and through the paper 

vf = y/'^ = K , (6) 

will be considered. We note that at the level of one-particle contributions, H, and as a conse- 
quence T and Vni, contain at the start all contributions entering in an elementary plaquette 
(unit cell for the system)^^. In these circumstances, for both Eqs.(4,5), it is important that 
r to be rigorously defined. This is because (i) we would like to represent different contribu- 
tions correctly in term of plaquette operators, and (ii) we must avoid multiple counting of 
different terms entering in the expression of H. For this reason, we mention that for a given 
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lattice site, taking into account nearest neighbor (NN) and NNN contributions as well (these 
will be present in an elementary plaquette), 8 hopping possibilities exist. Prom these, only 
4 are taken into account explicitly by J2n ^-iid, the remaining 4 contributions are introduced 
into the Hamiltonian by the H.c. operation. In these conditions, the defined 4 different r 
contributions entering in J2r 

X = axi , y = ax2 , x + y = a(xi + Xs) , y - x = a(x2 - Xi) , (7) 

where a is the lattice constant, and Xi, X2, are the versors of the Ox, Oy axis, respectively. 
The hopping and hybridization matrix elements generated by the contributions from Eq.(7) 
are represented for clarity in Fig. 1. Also for clarity, the explicit expressions of T and Vni 
from Eqs.(4,5) are presented in Appendix. A. 

In the case of concrete materials, the NNN one-particle contributions are small, and as a 
consequence are often neglected. Furthermore, it is important to know in rigorous terms if 
NNN contributions are introducing small corrections into the results or are able to provide 
qualitatively new effects. In the case of PAM, which is in a relatively early stage of its exact 
description, this issue must be also clarified. Because of this reason, during this paper, even 
if we start with NNN terms in H for technical reasons, we try to obtain exact ground-state 
solutions for PAM in the absence of NNN contributions as well. This task is also enhanced 
by the aim to extend the potential possibilities of the plaquette operator procedure we use. 
Starting from these motivations, we are reporting in this paper for the first time exact 
ground-states for PAM in 2D at finite and nonzero U, in the absence of NNN extension 
terms in H. 

In order to be able to obtain exact ground-states inD — 2 dimensions for the PAM Hamil- 
tonian presented in Eq.(l), we use a plaquette operator procedure which will be described in 
details in the following Section. Concerning the method itself, in our knowledge, it is used 
now for the second time (see also^^), and other methods in obtaining exact ground-states 
for PAM in D = 2 dimensions are not known at the moment. In principle, the procedure 
can be applied for other model Hamiltonians as well containing itinerant degrees of freedom. 
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The technique needs the transformation of if in a positive semidefinite expression based 
on plaqucttc operators. In 2D, the plaqucttc operators (as block operator units used for 
description) generate local, NN and NNN terms as well. This suggests that also the studied 
2D Hamiltonian must contain such type of contributions. If this would be the case, the ap- 
plication possibility of the procedure in 2D would be strongly limited to Hamiltonians that 
contain NNN extension terms as well. We demonstrate in this paper that this impression 
is not correct, and the procedure can be extended and apphed even in the absence of NNN 
contributions in H. We further mention that for Hamiltonians containing main long range 
terms (next to NNN or higher range contributions), the block unit used for description must 
be enlarged. 

III. PLAQUETTE OPERATORS USED FOR THE TRANSFORMATION OF THE 

HAMILTONIAN. 

Let us consider a 2D finite A^a — Nl x Nl square lattice, with lattice constant a. In 
order to identify the lattice sites, we are numbering them by i starting from the left-down 
corner, taking into account first the lowest row, and inside a row counting from left to right 
(we mention that for a vectorial position notation we are going to use i instead of i, when 
this is necessary). For example, in the simple case of = 4, we obtain the lattice site 
numbering presented in Fig. 2. As can be seen in this figure, we are denoting by Pi the 
elementary plaquettes. Using this notation, we start from the lowest elementary plaquette 
row, counting from left to right inside a row, and then going upward with the notation. 

Taking periodic boundary conditions into account in both directions, the number of 
plaquettes becomes equal to the number of lattice sites A^a- In this case, it is advantageous 
to denote every plaquette Pi by its down- left corner j, as pj. Concerning the notation of a 
plaquette through its down-left corner, for clarity we mention that for example, in Fig. 2., 
the plaquette P5 defined by the lattice sites (6, 7, 10, 11) becomes pe, or, the plaquette P7, 
defined by the lattice sites (9, 10, 13, 14) becomes pg, etc. 
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Let us now consider for pedagogical reasons, some c] fermionic operators creating parti- 
cles on lattice sites within the system. In general, the q operators can be labelled also by a 
supplementary a index which contains all relevant quantum numbers as well (in the case of 
PAM, a — {a, g), where a I denotes the spin, and g — d, f the type of particle). In this 
Section, being interested in the presentation of the method, we are neglecting the a index 
for simplicity. If the reader understands how the procedure works, the presented relations 
can be easily generalized for Cj^^ ^ well. 

Using the q operators, plaqucttc operators can be constructed by a linear combination 
of Cj acting on the corners of an elementary plaquette. We are denoting the coefficients of 
this linear combination by Up^i^c, where p denotes the plaquette, i labels a given corner of 
the plaquette p analyzed, and c denotes the type of operator considered (this becomes the 
a index when is used instead of Cj) , respectively. For example, in case of the plaquettes 
p{il) and p{il + 1) from Fig. 3A. we obtain 

Ail — (lp{ii),ii,cCil + Qp(il),a+l,cCil+l + Op(ii)jl,cCjl + Qp(il) jl+l,cCjl+l , 

Al+l = (lp(il+l),il+l,cCil+l + 0,p(il+l),il+2,cCil+2 + (ip(il+l),jl+l,cCjl+l + ap(il+l),jl+2,cCjl+2 • (8) 

Working with plaquettes, we must observe that all one-particle contributions of a given 
Hamiltonian can be obtained starting from plaquette operators. For example, let us consider 
the hopping matrix element connecting the nearest- neighbor lattice sites (i 1-1-1, j 1-1-1) from 
Fig. 3A, namely Tji+iji+i = (^ii+iji+icJi+iCji+i + H.c). This Hamiltonian contribution 
can be obtained for example, from the expression ^Ji^a -|- Indeed, we have 

ipa + ili+iia+i = Ta+iji+i + O , (9) 

where, the operator O concentrates all the other terms obtained from the left side of Eq.(9). 
The operator O will not be neglected in our considerations. It contains 30 terms that can be 
easily calculated from Eq.(8) (see also Appendix B.). The important aspect here, which must 
be keeped in mind, is that O do not contains contributions entering in Tji+i ji+i. Otherwise, 
the concrete expression of O is not important at the moment. The relation from Eq.(9) is 
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obtained since, Aj^An gives rise to (a*(.i) .i^i^^ap(a)ji+i,ccji+i9i+i + -H".c.), and 
creates the term {al^^^- i+i),ii+i,c'^p(*i+i),ii+i,cCii+iCji+i + H.c), respectively. Because the bond 
(il + 1, + is not present in other elementary plaquettes, even if we take into consideration 
all the plaquettes from the whole lattice in a sum of the form Z^j^J^j, the hopping matrix 
element t'ii^i ,ji_^^ becomes unambiguously expressed as 

^a+iji+i ~ ^p(ii),a+i,c^p(ii)ji+i,c + Op(a+i),a+i,c^p(ii+i)ji+i,c ■ (10) 

The obtained Eq.(lO) shows that the Hamiltonian parameters (at least the one-particle once 
in the present case) , can be expressed in term of plaquette operator parameters if we succeed 
to express the corresponding Hamiltonian terms into a sum of the form X^i ^l^i- 

Similarly, the next-nearest-neighbor hopping amplitude for the (il, jl + 1) hopping 
from the plaquette p{il) of Fig. 3A., contained in the Hamiltonian term Taji+i = 
(^iiji+i^iCji+i + H.c.) becomes 

= ^p(a),a,cS(ji)Ji+i,c ■ (11) 

In Eq.(ll) only the plaquette operator product ^Ji^a contributes, because, the NNN hop- 
ping described by Tji ji+i is contained only in the plaquette p{il). These examples illustrate 
that plaquette operators can be extremely useful in the study of different model Hamiltoni- 
ans H, since as seen from Eq.(9), different emerging contributions in H can be represented 
in diagonal, or positive semidefinite form via the operators Ai. As can be observed from 
Eqs.(10,ll), a such type of representation in term of plaquette operators, from the point 
of view of H parameters, simply means a parametrization in term of plaquette operator 
coefficients Up^i^c- For this to be possible, the plaquette operator products summed up over 
lattice sites of the form J2i A\Ai (i) must generate terms present in or (ii) must generate 
terms that are constants of motion (for example total number of particles, or lattice sites), 
or (iii) must generate terms that can be cancelled out if the (i) and (ii) conditions cannot be 
applied. We will return back to this problem after presenting the new plaquette operators 
defined in this paper (see after Eq.(17)), and the following Section exemplifies in detail a 
such type of transformation. 
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When the one-particle H parameters are not local (for example tl^^ji — ^a+iji+i — 
^ii+2,ii+2 ^ vertical nearest- neighbor hoppings), which means 

^a+iji+i = ' (12) 

the parameters j c of different plaquette operators are not independent. In the case of 
translational invariant Hamiltonians, we can chose for example translational invariant pla- 
quette operator parameters as illustrated by Fig. 3B. Denoting the sites inside a given 
plaquette starting from the down-left corner and counting anti-clockwise, the corners of the 
plaquette p{il) (and p{il + 1)) in Fig.SB., will be denoted by n (and m), respectively. Given 
by the considered translational invariance of plaquette operators, the plaquette operator 
parameters of the plaquettes p(il) and p{il + 1) with n — m — T equal indices will have the 
same value ar,c, t = 1,2,3,4. This property is extended as well to all plaquettes. In the 
examples contained in Fig. SB., the plaquette operators A^, Aji+i become in this case 

Ail — ai,cQl + Cl2,cCil+l + Cl4,cCjl + 03,cCjl+l , 

— ai^cQl+l + fl2,cCil+2 + fl4,cCjl+l + (l3,cCjl+2 ■ (13) 

From Eqs. (10,12,13), the unique NN hopping matrix element in y direction of c particles, 
based on Eq.(lO) becomes 

fy = a*2^^a3,c + a^^ai^c , (14) 

and, from Eqs. (11,13), the unique NNN hopping of the same particles along the main diag- 
onal of every elementary plaquette will be described by 

^l+y = «l,c«3,c ■ (15) 

Similarly, all one-particle Hamiltonian matrix elements can be expressed in term of plaquette 
operator parameters. When the so obtained equations (as Eqs. (14,15)) allow solutions for the 
plaquette operator parameters (this is possible usually in a restricted parameter space region 
Vh determined by the values of H parameters), the one-particle part of the Hamiltonian 
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can be expressed via (see Eq.(9)). Based on these relations and using for example 

properties related to positive semidefinite operators, the Hamiltonian of the system can be 
diagonalized exactly, at least for the ground-state, inside Vh- 

After testing this method in ID^^'^''' (using bonds instead of plaquettes), a such type of 
procedure has been recently used by us^^ in order to provide the first exact ground-state 
wave-functions for the periodic Anderson model in 2D in restricted regions of the parameter 
space. This has been done by choosing q = di^a-, fi,a for d, f electrons with fixed spin in PAM, 
and defining based on this choice, the Ai^„ spin-dependent plaquette operators containing 
spin-independent an,g parameters with g = d, f , n = 1,2,3,4 as follows (the example is 
taken for the plaquette p{il) of Fig. 3B.) 

= (il^dii^cr + Cll,ffil,a + a2,dC^a+l,(T + 02,//il+l,CT + 

0'3,ddji^l^o- + a^jfjl+l^fj + a^^dji^cr + 0-4,ffjl,a ■ (16) 

The obtained ground-state solutions based on Eq.(16) were connected to 3/4 fiUing^^, 
and are highly non-trivial states. One of them is a completely localized state, and the 
second one is itinerant, with the momentum distribution function for the half- filled upper 
diagonalized band as shown in Fig. 4, presenting a clear evidence of (exactly deduced) non- 
Fermi liquid behavior in normal phase and D = 2 spatial dimensions. This shows that the 
procedure detects ground-states which are far to be trivial. However, the inconvenience of 
the plaquette operator from Eq.(16) is that via J2i,cr M,a'^i,a it creates NNN terms, these 
must be present in H as well, so the deduced ground-states, and the procedure itself, seem 
to be related to the presence of NNN extensions in the Hamiltonian. We present below how 
this inconvenience can be removed. 

For this reason, we must observe, that the choice of the operators q in Eq.(13) and the 
form of the plaquette operator itself is not fixed a priori. This means that the possibility 
presented in Eq.(16) for the plaquette operators is not unique, even if we are interested 
in the study of a fixed model (as PAM in the present case). As a consequence, we can 
chose other possible forms for the plaquette operators, and using them, we can deduce other 
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ground-states in other regions of the T — phase diagram of the model. To exemphfy 
this statement, in the present paper we define for the decomposition of the studied PAM 
Hamiltonian in translational invariant case, a completely new type of plaquette operators 
Ai and Bi. Each of these has different plaquette operator parameters an,g,a and bn,g,a, 
n — 1,2,3,4, g — d,f, a ==t) i- Furthermore, both plaquette operators Ai and Bi are 
containing both spin components with different numerical prefactors, i.e. an,g,a, bn,g,a are 
considered independent, and a dependent. Exemplifying the new form for the case of the 
plaquette p{il) of Fig. 3B., where i denotes the vectorial position of the site il, the new Ai 
operator is defined as 

Ail = ai,d,|iiii,T + 02,d,|iiii+i,| + 03,^,1^^1+1, | + a4,d,|(iji,| 
+ ai,d,idii,i + a2,d,idii+i,i + a3^d,idji+i,i + «4,da^ji,i 
+ ai,/,T/ii,T + 02,/,T/ii+i,T + a3,/,T/ii+i,T + «4,/,T/ii,T 

+ aij,ifii,i + a2,f,ifii+i,i + asj,ifji+i,i + aij^Jji^i . (17) 

Similar expression is used for the Bi operator as well for the same plaquette p{il), in which, 
the plaquette operator parameters are considered &n,fl,o-) instead of an,g,n- Note the plaquette 
independent values of the an,g,<T and &n,g,<T parameters, which, as explained in this Section, 
is given by the translational invariance of the considered system. 

Comparing Eq.(16) and Eq.(17), we realize that the Aj,^. plaquette operators for both 
(T =t, J, values have 8 independent a„,g parameters, while in the present case, for both 
ylj, Bi operators, the number of independent plaquette operator parameters is 32. This 
enlargement of the number of parameters give us the possibility to demonstrate that the 
described procedure is able to detect also ground-states whose presence do not require the 
NNN terms in H of the system, even if the A\Ai products are providing such type of terms 
at the start. The key feature for this to work is the presence of two plaquette operators Ai 
and Bi containing different spin-dependent coefficients. Indeed, in this case, by B\Bi we 
can cancel out not only the || terms created by Y^iA\Ai which arc not present in H (these 
would represent for example hopping terms containing spin- flip), but also the NNN terms 
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generated by AjAi. Because of this reason becomes to be possible to obtain the expression 
of a Hamiltonian not containing NNN contributions in term of plaqucttc operator products 
which create such type of elements. The concrete transformation of the Hamiltonian is 
presented in the following Section. 

IV. THE HAMILTONIAN WRITTEN IN TERM OF PLAQUETTE OPERATORS. 

Comparing the expression of the Hamiltonian presented in the previous Section together 
with the exphcitations contained in Appendices A. and B., Eqs.( 3,A1,A2,B1), we reahze 
that the following relation holds 

Na Na 
i i 

4 Na 4 Na 

+ E[E(l«n,.,.r + l&nA.r)] E4A. + E[E(l«n,/,.r + l&n,/,.n] E/i!X , (18) 
f n=l i c n=l i 

if, the hopping and hybridization matrix elements are related to the parameters of the 
plaquette operators Ai and Bi via 

F{t',,Vr;an,g,a,bn,g,a)=0, (19) 

where the non-linear system of equations from Eq.(19) is presented explicitly in the Ap- 
pendix. C, and g — d,f. These equations arise as Eqs.(14,15) in Section I. The system of 
equations Eq.(19) must be considered as containing known H parameters (t^,K), and un- 
known plaquette operator numerical prefactors {an,g,a, i>n,g,a)- In fact, a simple (but lengthy) 
algebraic calculation shows that Eq.(18) exactly holds if the relations between the param- 
eters of H and the numerical prefactors of the plaquette operators, presented explicitly in 
Appendix C, are satisfied. The number of equations contained in Eq.(19) is 70, and the 32 
unknown complex plaquette operator parameters provides 64 unknown variables (the real 
and imaginary parts). These are entering in Eq.(19) in a nonlinear, but complex-algebraic 
manner. Since the number of equations is higher that the number of unknown variables, 
solutions will be allowed only if some inter-dependences (fixed by Eq.(19)) will be present 
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between the H parameters. These relations contribute in the definition of Vh (see also the 
observations below Eq.( 26)). 

We underline that since the structure of plaquette operators used in this paper (see 
Eq.(17)) is completely different from the structure of the plaquette operators from Eq.(16) 
(in that case, instead of Eq.(19), we have 17 equations presented in Eq.(9) of Ref.^^, con- 
taining 16 unknown variables), the problem set up here, from mathematical point of view, 
is completely different from that analysed in our previous work. 

We also note that as can be seen from Appendix B, X^i ^i^i introduces (t, i) like terms 
as well, which are missing from the Hamiltonian. Because of this reason we need a second 
plaquette operator product BjBi, whose role is exactly to cancel out these supplementary 
terms not present in H, Eq.(l). Furthermore, the presence of Y.\BlBi allows also to cancel 
out the NNN terms created by X^i^i^i. Via Eq.(18), this give as the possibility to express 
H through Ei(4^i + B\B-;) even in the absence of NNN terms in the Hamiltonian, and to 
obtain ground-state wave-functions in this case as well. 

Using now Eqs.(2,3), we have 

lJ + Ef = UP' + {Ef + U)Y: flka - UNa , (20) 

i,<T 

where, the positive semidefinite operator P' has been defined in Sec. II. Adding Eq.(20) to 
Eq.(18) and using for the plaquette operators the anti-commutation property presented in 
Eq.(B2), we find 

H = + BiBl) + UP' - Na{U + K'l + Kf + Kf + K{) 

i 

+ N^K^ + NfKf + iV/i^f + iV/i^f + {Ef + U) (iV/ + TVf) , (21) 

where, the introduced constants are defined by — Y^n=i{\^n,g,(j\'^ + l^n.^.o-P)) and the 
particle number operators by 7V| = Si^^' with g — d, f. Imposing the relations 

K^^Kf^K, K\^kI^K^, Ef + U^K-Rf, (22) 

the expression of H from Eq.(21) becomes 

14 



H = Y.iAA\ + B;BI) + UP' + KN- Nk{AK - 2Ef - U) . 



(23) 



Since we are working at fixed number of particles N , from Eq.(23) we obtain 



P + E, 



(24) 



where P for [/ > is a positive semidefinite operator defined by 



P = ^(iii[ + SiS|) + C/P', 



(25) 



and the number Eg is given by 



Eg^KN - Nj^{4K - 2Ef - U) . 



(26) 



The transformation of Eq.(l) into Eq.(24) is possible only if the system of equations 
Eqs.(19,22) allows solutions for the plaquette operator parameters. The presence of these 
solutions will be possible only on restricted domains Vh of the parameter space of the prob- 
lem given by the inter-dependences between the H parameters mentioned below Eq.(19). As 
a consequence, the solutions that will be presented below are valid only in this Vh region. 

V. EXACT GROUND-STATE WAVE-FUNCTION SOLUTIONS. 

In this Section we are presenting first the derivation of the exact ground-states, then we 
discuss the possible solutions for the plaquette operator parameters, and finally, we analyse 
in extreme details the solution obtained for zero NNN contributions. 



Starting from Eq.(24), taking into account that P is a positive semidefinite operator, we 
realize that the ground-state of H = P + Eg is the wave function |^c,), for which P|^'g) = 0. 



AlAl = , BlBl = , AlAl = -AjAl , BlBj = -BjBl , AIb} = -BjAl , (27) 



A. The derivation of the exact ground-states. 



To find l^g), we have to keep in mind Eq.(25) which defines P. Given by 
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we observe that the plaquette operator part of Eq.(25) apphed to Hi^i-^/ gives zero. Fur- 
thermore, since P' requires for its zero (and minimum) eigenvalue at least one /-electron 
on every lattice site, we add to the ground-state the contribution = ni(A*i,T/i^t + l^hlfli)^ 
where /li^a ^-re arbitrary coefficients. As a consequence, the ground-state with the property 
Pl'^g) = becomes 

1^.) = U[MBl{fii,^fl+fiiJi)m , (28) 

i 

where, |0) is the bare vacuum with no fermions present. The product in Eq.(28) must be 
taken over all lattice sites. Because of this reason, the product of the creation operators 
in Eq.(28) introduces N — 3Na particles within the system, so the deduced ground-state 
wave-function corresponds to 3/4 fiUing. All degeneration possibilities of the ground-state 
are contained in Eq.(28), since the wave function with the property Pl'if) = at 3/4 filling 
always can be written in the presented form. We underline however that PAM contains 
two hybridized bands, and 3/4 filling for a two-band system means in fact half filled upper 
hybridized band (the lower band being completely filled up). 

The wave- vector \^g) represents the ground-state of the starting Hamiltonian, only if 
Eq.(l) can be transformed in Eq.(24). This is possible only if we are situated inside the 
region Vh of the parameter space, i.e. the system of equations Eq.(19) detailed in Appendix 
C. allows solutions for the plaquette operator parameters, in conditions in which also the 
constrains from Eqs.(6,22) hold. In the remaining part of the paper we will concentrate on 
these possible solutions. 

We underline, that presented in Eq.(28) describes rigorously only the U > case, 
since the presence of the Fj^ operator into the ground-state is required only by the non-zero 
U value. As a consequence, the ground-state at U = cannot be expressed in the form 
presented in Eq.(28). We emphasize that the differences between Eq.(28) and the ground- 
states deduced previously^^ are present because instead of Y[i,a ^l,a obtained in the old case 
with Ai cr defined by Eq.(16), we now have in the ground-state wave function Hi^i-Bp 

Before going further, we mention that the physical properties of the ground-state wave- 
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function written mathematically in Eq.(28) strongly depend on the nature of the concrete 
solutions provided by Eq.(19). 



B. Solutions for the plaquette operator parameters. 



The solutions for the plaquette operator parameters which lead to the ground-state 
\^g) must be obtained solving together Eqs.(6,22,Cl). These taken together represent 74 
nonlinear complex-algebraic coupled equations, so a relatively difficult mathematical task. 

A study of the next-nearest neighbor contributions entering in Eq.(Cl) shows that the 
solutions exist only if the following inter-dependences are present between the plaquette 
operator parameters 



where, x^, Vu are complex, finite, nonzero, otherwise arbitrary parameters defined by the 
ratios presented in Eq.(29). Using Eq.(29), the studied system of equations can be completely 
transcribed for the &n,g,(T unknown variables with n — 1,2,3,4; g — f, d; a —'\, \ (the 0^,9,0- 
parameters being given through hn^g^a via Eq.(29)). Since the so obtained equations for the 
bn,g,c variables are representing the starting point of the description of physical properties 
provided by the deduced ground-states, they are presented in Eqs.(Dl,D2) of Appendix D. 
Starting from this moment, we must solve the system of equations presented in Appendix 



We have found for the system of equations Eqs.(Dl,D2) several mathematical solutions, 
which will be briefiy presented below. 

a) Taking x^ — — y, xi — y^ — —1/y*, we find the first class of solutions. The 
interesting aspect of this case is that the 20 equations contained in Eq.(Dl) are automatically 

satisfied, and we must concentrate only on equations presented in Eq.(D2). This last system 





(29) 



D. 
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provides a solution for bn,g,i — ybn,g,'i, 9 — d, f, which however do not has new aspects in 
comparison to the solutions we find in Ref.(^^). 

b) As can be seen, in order to obtain new solutions, we must take ^ y^, xi ^ yi into 
account. The first attempt that can be made, is to consider xi — xi — x, yi — yi — y, and 
X — y. This solution presents the interesting property that reduces the system to ID case. 
This means that the solution emerges only for = = = and t^^^ = tyj^^^ = Vy±x = 0. 
New aspects related to PAM in comparison with those reported in Refs.(^^'^^) are not present. 
This case merits however attention in the future, since it allows to study at the level of exact 
ground-states (taken in the form of Eq.(28)) the modification of the ID properties to 2D 
characteristics by taking into account small and smooth deviations from the x — y condition. 

c) The third solution that we have found was deduced in = = x, yi = yi = y, 
and X ^ y case. This solution will be presented here in details, since presents a 2D ground- 
state that emerges for zero next-nearest-neighbor H contributions. A such type of exact 
solution for PAM is completely new, because it cannot be obtained by the decomposition 
used previously^^. 

d) We have studied also the general x^ 7^ x^^, y^', y^ 7^ y-a, ^a'l case as well, obtaining 
only localized solution which require the presence of next-nearest neighbor H terms as well. 

C. Detailed analysis of the solution obtained in the absence of next-nearest neighbor 

Hamiltonian terms. 

Herewith, we analyze in detail the solution c) described above requiring x ^ y. This 

emerges at 

^y±x ~ ^y±x = ^y±x = ^y±x = ) (30) 

SO it describes a ground-state wave function for PAM not containing in its Hamiltonian 
NNN extension terms. A such type of exact ground-state in 2D at finite nonzero value of 
the interaction is presented for the first time in this paper. 
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Solving for the plaquette operator parameters the system of equations Eqs.(Dl,D2) we 
have found 

* yx*r2 y* 1 

ai,d,T Pd, = Pd ' «2,d,T = —Pd , a2,d,l = Pd , 

Zl Zl Ti 

_^2|/* _ _ 1 _ 1 

Os.d.T — —Pd : 0'3,d,l — —Pd : 0'4,d,'[ — jPd , 0,4^4,1 — Pd ) 

Zl yTi Zl 



yx*T* y* 1 
«i,/,T = ^*Pf ' = ^^P/ ' «2,/,T = -^Pf , 0.2 J,i = -—Pf , 

Zl Zl Ti 

T2y* 1 1 

^1 yT^i Zl 

Ol,d,T - ' "l,d,l - —Pd > ^^2,0!,^ - —pd , 02,d,l " "TTTPd , 

-Zl -Zl 'ly 

xT2y* , 1 y 

03,d,T = ~Pd , 03,(i,i = , 04,d,T " "^'ci ' ^"^Al — — Pd ) 

Zl Ti Zl 

KfA = Pf > Kf,i = —Pf > 02,/,T = -jP/ , 02,/a = -rrrP/ > 

Tiy* 

:^Pf > ^3,/,i = , 64,/,! = —Pf , b4,f,i = —Pf . (31) 



" n Zl 

The conditions imposed for the parameters entering in Eq.(31) are xy* ^ —1, x ^ y, ti^ T2, 
T1T2 = real, pdP*f = real. Together with Eq.(31), the nonzero H parameters become 

4 = -Ri\pd\', ti = -R2\pd\\ ti = -Ri\pf\\ tf = -R2\pf\\ 

V,^-Rip*^f, Vy^-R2P*dPf, Vo^-Rsp*dPf, U + Ef^K-Kf, 

K = Rs\pd\\ Ki = R^\pf\\ (32) 

The Rni n = 1,2,?) factors present in this relation are given by 

zl Ti yrt \ziY 

fi3 = (l + |.f)(l + ^) + if^(l + ^). (33) 

We further mention, that the obtained solution, for y = zi = ti and x = T2 reduces to the 
isotropic case, where — — , g — d, f , and — Vy = V. 

The ground-state wave function from Eq.(28) in the case of the solution from Eq.(31) 
reduces to the simple form 
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i 

The result presented in Eq.(34) is obtained because A = studied case, is 

unable to introduce three particles on the same lattice site. Since Pi introduces N — 2Nj^ 
electrons in the system, being unable to put three electrons on a given site, an uniform 
particle distribution will be obtained with two electrons per site, which generates the product 
Oibd^^i,! + P}fl]){P*ddli +P*ffli) in l^gi^oc)) in Eq.(34). The added term contained in 
(see Eq.(28)) introduces one more / electron on each site, and as a consequence, do not 
modify the created uniform particle distribution within the system, and Eq.(34) arise. This 
wave-function has a well defined norm 

i 

and as mentioned above, coherently maintains three particles on every lattice site 

h,\^b,(loc)) = E(4.4. + flA'rmAloc)) = S\^,(loc)) . (36) 

Denoting by (...) = g{loc)\...\^ g{loc)) / g{loc)\^ g{loc)) the ground-state expectation val- 
ues, we obtain long-range density-density correlations within the system 

/^(niUj^i) = 3. (37) 

Furthermore, it can be observed that (/oc)) prohibits in the same time the hopping and 
non-local hybridization between all site pairs 

(f,) = {ff) = {V^i) = , (38) 

since, for g, g' = f, d, we have {gj^f^g'jy) = 0, for all j 7^ i. As a consequence, the ground-state 
l^p(Zoc)) clearly represents a completely localized state. 

The remaining non-zero ground-state expectation values of different H terms are given 

by 

{Efl^ 2\p£+\p£ 0) _ \pf\' (Vo)_ K\pf\' + Kf\p,\' 

Na ^ Ip.P+Ip/P ' \Pd\' + W Nj, \Pd\' + \Pf\' ■ ^ ' 
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Because of C/ > 0, and as seen from Eqs.(32,33), K, > 0, the non-zero on-site hybridiza- 
tion, couphng the two bands, decreases the energy of the system ((t^o) < 0, (£"/) > 0, {U) > 
0). The ground-state energy becomes Eg = (Vq + Ef + U) . 
In average, the number of /-electrons per site becomes 



(E/U<^> = i + Y^. (40) 

where, x = Pf/pd and we have x"^ = It'^ /f^l in the isotropic case, and x^ = \t{^/ti\, a = x,y 
in general, for the considered solution. Since in concrete physical situations x « 1, the 
number of / electrons per site is close to one, but not exactly one in the ground-state. 
Excepting the small number of sites with double /-electron occupancy, the local /-moments 
are not compensated. In fact, defining mf — nf-^—nf^^, with g — d, f, we have (ifif) — x{m,(), 
and (mf) -|- (mf) — (|A<i,-f|^ — |A*i,j,n/(|A''i,Tp + l/^iaP)' where as presented before, //i^o- are 
arbitrary. 

Concentrating on the magnetic properties of the ground-state, the total spin of the 
system can be standardly expressed via — l/2J2iidl^idi^\ + fl\fi,] ~ ~ /u/i,i)' 

= J:i{di^di,i + fl^ki), S- = {S+y, and = Si + 1/2(5+5" + S-S+). Calculating the 
ground-state expectation values, we find 

/C2^ ^ SA^A ^ (|/Xi,tP - |/^Up)(|/^j,Tp - l/^jgH + 2(//i,|/i?^/ij,t/X?^ + C.C.) 

^ ^ ^ k 4(K,P + K,P)(|/.j,,P + K,P) 



Taking now two extremum {fJ^i^a} distributions, a) for /^i^i = /j,^, Hi^i = 0, we find 
(Si) = {Na/2Y, and {S'^) = {Na/2){Na/2 + 1). This situation corresponds to maximum 
total spin in the system, with average total spin absolute value per site {{S'^) / Nj)^/'^ — 



1/2(1/2 -|- I/JVa), which is of order 1/2 for large A'a. b) Dividing however the square lat- 
tice into two equal sub-lattices with //i | = //i ^ = in one sub-lattice, and /Xi^^ = 0, /Xi^j^ = ^ 
in the other one, we obtain {Sl) — 0, \f{S^)/Nl — 1/(-\/27Va), i.e. in the thermodynamic 
limit, the total spin in absolute value per site is zero in this case. As can be observed, 
the degeneration of the ground-state physically is given by the fact that all possible total 
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spin values are contributing in its construction. As a consequence, the ground-state behaves 
paramagnetically. We must further observe, that not all different {/ii,cr} sets provide linearly 
independent ground-state wave-function contributions. For example, choosing for all i the 
values /li^i — III, /li^i — 0, or /li^^ — 0, /li^i — /i^, or /li^^ — /i-^, /li^^ — /ii, we recover the same 
ground-state with maximum value of {S'^). As a consequence, in order to find orthogonal 
wave-functions that belong to the ground-state, the /^i^o- coefficients cannot be chosen com- 
pletely random and independent. Also the normalization to unity of |^p(Zoc)) represents a 
constraint to the value of these coefficients. Neglecting the trivial multiplicity obtained from 
the spatial orientation of the total spin S, the degree of the degeneration of the ground-state 
is NJ2. 

The spin-spin correlation functions can be also calculated, and for i 7^ j we obtain in the 
case of a fixed {lJ>i,a} set 





2 






/'j.lP 


2 _ 


/'j,J 




-2(// 


lift 




4 




( 


l/^i,T 







(42) 



Since, as shown before, the Hi^a coefficients are not completely independent, the spin-spin 
correlations given by Eq.(42) are quasi-random. Resembling behavior is experimentally seen 
in heavy- fermion cases^^. 

The phase diagram region where the solution occurs is presented in Fig. 5. for the 
isotropic case. The general aspect of this region remains the same in the anisotropic case as 
well. It represents a surface in the parameter space which extends from the low U region 
up to the high U region as well. This region is completely different from that obtained in 
Ref.(^^) which emerges for nonzero values of next- nearest- neighbor hopping and non-local 
hybridizations. 

The non-local nearest-neighbor hybridization matrix element in the isotropic case is 
related to hopping matrix elements by (y/i*^)^ = /t^. In the anisotropic case this relation 
becomes (K/^a)^ = ti/t^, a — x,y. Modifying the values of hopping or/and hybridization 
matrix elements, we can leave Vh, destroying the ground-state character of |^'g(/oc)). This 
process can be tuned by pressure which strongly influences the U, Vr parameters (see for 

22 



example Ref.(^^)). Since the reduction of Eq.(28) into the completely localised |^g(Zoc)) 
from Eq.(34) it is itself based on a dehcate balance between H parameters (contained in 
Eqs.(32)), the loss of the localization character of particles in principle can be easily achieved. 
This localization-delocalization transition represents in fact a MIT transition provided by 
PAM. Since the exactly deduced ground-state energy do not contains the exponential term 
characteristic for a Kondo type behavior (see for example the discussion presented in Ref.(^)), 
a such type of MIT transition cannot be connected to Kondo physics. Instead, the MIT 
transition connected to the destruction of the localized \^ g{loc)) ground-state is related to 
the break-down of the long-range density-density correlations. 



VI. SUMMARY AND CONCLUSIONS 

We are presenting for the first time exact ground-states for the periodic Anderson model 
(PAM) at finite U mD = 2 dimensions in the case in which the Hamiltonian does not contain 
next to nearest neighbor (NNN) extension terms. For this reason, and based on this frame 

(i) The used plaquette operator procedure is presented in detail and it is underlined that its 
applicability is not connected to the presence of NNN extension terms in the Hamiltonian. 
We underline that this is the only one procedure known at the moment, which is able to 
provide exact ground-states for PAM va. D — 2 dimensions and finite U interaction values. 

(ii) A new plaquette operator has been introduced for the study of the PAM Hamiltonian. 
The plaquette operator contains contributions coming from all spin components, possesses 
spin dependent numerical prefactors, and allows the detection of ground-states even in the 
absence of NNN extension terms in the Hamiltonian in restricted regions of the parameter 
space, (iii) The physical properties of the deduced ground-state have been analysed in detail. 
All relevant ground-state expectation values and correlation functions have been deduced 
for this reason, (iv) The implications of the deduced results relating the metal-insulator 
transition in the frame of PAM have been analysed. It has been pointed out that the 
lost of the locahzation character in the studied case is connected to the break-down of the 
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long-range density-density correlations rather than Kondo physics. 

The obtained new exact ground-state emerges at 3/4 filling, not requires the presence 
of NNN extension terms in the Hamiltonian, it is paramagnetic, and presents quasi-random 
spin-spin correlations. It represents a fully quantum-mechanical state (in the sense that it is 
far to be quasi-classical), and it is build up through superposition effects. The ground-state 
wave-function coherently controls the occupation number on all lattice sites, introducing 
in this manner long-range density-density correlations within the system, and prohibiting 
in the same time the hopping and non-local hybridizations. The local / moments are not 
compensated and the / electron occupation number per site in average is close to, but not 
exactly one. 

Concerning the question of the physical relevance, we would like to mention that in 
general terms, even a solution detected in a restricted parameter-space region which behaves 
completely repulsively in the renormalization group language could have significant physical 
implications^^. In the present case, besides presenting open roots toward the deduction 
possibilities of exact ground-states in D = 2 dimensions for strongly correlated systems, 
the presented results provide exact theoretical data which can be used in the process of 
understanding and description of the metal-insulator transition in the frame of PAM. 
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APPENDIX A: THE EXPLICIT EXPRESSION OF THE NON-LOCAL ONE 
PARTICLE CONTRIBUTIONS IN THE HAMILTONIAN. 



The kinetic energy term for d electrons has the expUcit form 
Td = E[(^'4<^^"i+x,. + H.c.) + {tidlA+y,, + H.c.) 

i,(7 

+ itUyd\,A+{^+y),a + H.C.) + + H.C.)] . (Al) 

The kinetic energy term for / electrons can be simply obtained from Eq.( Al) interchanging 
the d index with the / index, and d by /. 

The explicit expression of the non-local hybridization becomes 

Vnl = E[(^/4. W + H.C.) + {Vj^flA+^,a + He.) 

+ {VfdlJ,+y,, + H.C.) + {Vj'flA+y,, + H.C.) 

+ {VSydlMi^+y),a + H.c.) + {VityflAH^-^yU + H.c.) 

+ Ji+(y-x),. + H.c.) + {yf^JlAHy-^),a + H.c.)] . (A2) 



APPENDIX B: THE PLAQUETTE OPERATOR CONTRIBUTIONS SUMMED 

UP OVER THE LATTICE SITES. 

The expression ^^Ai summed up over the whole lattice considered with periodic boundary 
conditions in both directions is presented below in condensed form (^f, g' = d, f ; g , g' = d, f). 

Na 

i a,(T' g,g' i 

[^iV^i+x,(T'(«t,g,cT«2,p',a' + <g,aa3,g',a') + H.C.] + [gij[+yy{a\^g,^ai^g>,a> + a\g^^a^,g> ,a>) + H.c] + 

[^iV^i+(x+y),<T'(«t,s,a«3,9',<T') + H.c] + [^iV^i'+(y-x),a'(«2,5,a«4,5',a') + H.c] + 

[^U>'(E <,g,a<^n,g',a') + H .C.][l - Ug^g^d^A) ■ (Bl) 
n=l ^ 

Furthermore, the following property is satisfied 
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A\A, + A,A\ = Y.{\an4A? + Kai? + Kja? + Kja?) ■ (B2) 



n=l 



For the plaquette operators the Eqs.(Bl,B2) hold as well by changing the coefficients 
an,g,<j to 6n,fl,a, where g = dj. 

APPENDIX C: THE NONLINEAR SYSTEM OF EQUATIONS. 

The explicit expression of the system of equations Eq.(19) containing 70 equalities is 
presented below in condensed form. The used abreviations are g,g' = d,f, a =T,i, and 
(g^gi) represents (dj) or 

= ^\,g,(T^2,g,a + (H,g,a^2,,g,a + b\,g,<T^2,g,a + bl,g,ah,g,a , 

= ^l,g,(T^4,g,(T + 0'2,g,a'^3,g,a + ^l,g,cr^4,g,<T + ^2,g,cr^3,g,<T , 

= ai,5,(T«2,g' -a + al g^^as^g'-a + bl g^^b2,g'-a + hl g^^h2„g>-a , 

= «l,g,o-«4,g' -0- + a^^g ^^U^^g' -a + h\ g j}4:,g' -a + ^'2,g,cr^3,g' -o" , 

= Ql,g,o-^3,5,(T + W,g,a^3,g,a , —t^-x = '^2,g,a'^4,g,a + ^2,5,o-^4,p,(T , 

= ^l,g,a(^3,g',-a + &t,g,a^3,s' -f^ ' = «2,g,<T«4,s' -a + ^'2,fl,a^4,g' -a , 

-V/'^^ = Oi^g o.a2,gi,(T + a4_g_^a3,gi,o- + bl,g,<j^2,gi,a + ^'4,fl,a^3,fli,o- , 

-V^jf'^^ = «l,g,(T«4,gi,(T + 02,3,(^03,91,(7 + K,g,c^^,gi,(T + ^2,g,(T^3,gi,(T , 

"l^x+j/ — '^l,g,a^3,gi,a + '^l,g,a'^3,gi,a , ~^y-x — '^2,g,a^'^,gi,(T + ^2,g,a'^4:,gi,a , 
4 4 

n=l n=l 
4 4 

n=l n=l 

APPENDIX D: THE EQUATIONS FOR THE PLAQUETTE OPERATOR 

PARAMETERS 



After using Eq.(29), the remaining equations for the plaquette operator parameters are 
presented in detail in this Appendix. These equations can be divided in two parts: a 
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homogeneous part (Eq.(Dl)), and a non-homogeneous one (Eq.(D2)), see below. These two 
system of equations are presented here in detail in condensed form. Written explicitely, 
Eq.(Dl) (Eq.(D2)) contains 20 (41) different equations, respectively. 

The homogeneous part of the equations is as follows {a =t, i, g — d, /, g' — d, f) 

(-^ + m,,,M,,',-a + (— ^ + 1)&:,,,.&4,,',-. = , 

+ m,,,ah,,',-a + {-^ + = , (Dl) 

The non-homogeneous part of the equations is presented below. The abbreviations used 
here are {g,gi) — {d, f) or {f,d), and the presence of a means that two equations are 
simultaneously present with cr =t and cr =J,. For a single g index we have g = d, f. 



-ti = {X^y; + 1)6I,,,<,62,,,. + + m,g,ah9,cr , 

X—cy—cr 

y—a •I'—a 

X—ay—a 

y—a -^—(j 



(^ + l)l&3,.,.P + (^ + 1)164. 



X—a y—a 

-Vx+y = i~~ + ^)K,g,ah,gi,a , -Vy-x = (~~~ + ^W2,g,ah,gua , 
X—a y—a 

Vf^Vf^^Va, a^x,y,x + y,y-x; U + Ef ^ K - , K ^ . (D2) 



27 



FIGURES 




FIG. 1. Hopping and hybridization matrix elements indicated by dashed lines in the elementary 
plaquette {il,i2,i3,i4). 
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a 



PI 



P2 



P3 



a 

FIG. 2. A 4 X 4 square lattice in 2D covered by elementary plaquettes denoted by PI, 
I = 1,2, ...9. The numbers paced in down- left corner of every plaquette denote the lattice sites, 
and a is the lattice constant. 
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jl 



jl+2 



B) 




il il+1 il+2 

FIG. 3. Two neighboring plaquettes p{il) and p{il + 1) in the square lattice (A), and the 
notation of lattice sites inside the plaquettes by n and m respectively (B). 
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n(k) ' 
1 



0.5 



FIG. 4. Momentum distribution function n(k) for the itinerant solution of Ref.^^ in the upper 
diagonalized, half filled band, along the whole first Brillouin zone. As can be seen, neregularities of 
any kind in n{k) and its derivatives of any order are missing, signalling non-Fermi liquid behavior 
in normal phase and 2D. 
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3 




FIG. 5. Phase diagram region where the locahzed solution occurs in the absence of 
next-nearest neighbor terms and isotropic case. The presented surface extends up to infinity for 

{U + Ef)/\t'^\ ^ oo. 
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